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Abstract 

We examine the weak-field approximation of locally Galilean in- 
variant gravitational theories with general covariance in a (4+1)- 
dimensional Galilean framework. The additional degrees of freedom 
allow us to obtain Poisson, diffusion, and Schrodinger equations for the 
fluctuation field. An advantage of this approach over the usual (3 + 1)- 
dimensional General Relativity is that it allows us to choose an ansatz 
for the fluctuation field that can accommodate the field equations of 
the Lagrangian approach to MOdified Newtonian Dynamics (MOND) 
known as AQUAdratic Lagrangian (AQUAL). We investigate a wave 
solution for the Schrodinger equations. 
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1 Introduction to the Galilean framework 



It is well known that the Galilei group describes non-relativistic, or low- 
energy, phenomena pQ. The study of non-relativistic limits of various rela- 
tivistic phenomena is not always a straightforward procedure; for instance, 
the straightforward low-velocity approximation of electromagnetism may 
lead to one limit, whereas there exist in fact two limits of electromagnetism 
[2J Another example of such a subtlety is the concept of spin; although 
historically explained within the context of the relativistic Dirac equation, 
spin can also be understood in terms of the Galilean theory, so that it is not 
a purely relativistic notion [J]. 

The Galilean symmetry is, in some respects, more intricate than the 
Poincare symmetry, which underlies relativistic physics [JJ. For instance, 
Galilean systems are not usually described in terms of tensors, because of the 
absence of a metric. This problem was circumvented, for flat manifolds, with 
the formalism of Galilean covariance, originally examined by various authors 
[5]- [7]. This formalism consists in describing Galilean theories in a tenso- 
rial form by means of a light-cone metric defined on a (4 + l)-dimensional 
Minkowski manifold. As mentioned in these early papers, the reduction from 
this extended Minkowski space to a (3 + l)-dimensional spacetime can lead 
to the Galilean invariant theory as well as the Minkowski theory, thereby 
providing a unified framework to treat both the Galilei and the Lorentz 
kinematics in a simple and elegant manner [6]- [7]. 

In this paper, we extend the formalism to the Riemannian geometry 
underlying gravitational models, following the lines of Ref. [8]. One motiva- 
tion for this work is the fact that General Relativity (GR) does not provide 
satisfactory explanations to some gravitational phenomena involving non- 
relativistic objects, e.g. the dark matter and dark energy problems [9]. The 
presence of dark matter is required in order to explain, for instance, the 
rotation curves of spiral galaxies, and dark energy appears in the context 
of cosmology in order to describe the accelerated expansion of the universe. 
Although quite diverse, these issues present a common feature: they occur 
in the low-energy regime. This indicates that local Galilean invariance as 
the underlying symmetry in the theory of gravitation could shed new light 
on these questions. The MOdified Newtonian Dynamics (MOND) approach 
[10] -[13] addresses these issues either by modifying the law of inertia or by 
altering the Newtonian theory of gravity phenomenologically. Our perspec- 
tive is aligned with the latter approach; we will show that it accommodates 
the AQUAL fundamental equation in a rather natural way, with the extra 
dimension being related to the mass density. 
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The Galilean formalism, described in Refs. [5]- [TJ, consists in build- 
ing Lorentz-covariant action functionals defined on a (4 + f )-dimensional 
Minkowski manifold. This manifold is described in terms of the Galilean 
five-vector, (x, x 4 ,x 5 ), which transforms under Galilei boosts, with relative 
velocity v, as 



4 

x — vr, 



x' 4 = x 4 , (1) 
x' 5 = X 5 — V • x + iv 2 x 4 . 

Thus, general Galilean transformations (including rotations, Galilei boosts 
and translations) have a form similar to the usual Poincare transformations, 

x ,fl = A^x u + a" , fi,i/ = l,...,5. (2) 

The Galilei algebra of a (3 + l)-dimensional space-time is an 11-dimensional 
subalgebra of the 15-dimensional Lie algebra isomorphic to the Poincare 
algebra in (4 + l)-dimensional spacetime. 

The tensor methods utilized in the Galilean covariance formalism are 
based on the scalar product, A ■ B = A • B — A^B^ — A5B4, which is invariant 
under the transformations in Eq. (pQ). This suggests to employ, on a locally 
flat manifold, the following Galilean metric: 



3x3 














-1 





-1 






= -1 . (3) 

V -10/ 

We call Galilean manifold this 5-dimensional flat manifold which, although 
equivalent to a Minkowski manifold in (4 + 1) dimensions, reduces to the 
Galilean space-time via appropriate ansatze. The five-momentum, pn = 
idy = (1V ,idt,ids) = (p, —,mc\, where c is a parameter with the dimen- 
sions of velocity, and S is the extra coordinate, suggests that the additional 
coordinate, x 5 = — , may be seen as the conjugate of mc, where m is the 
mass. 

In this paper, we utilize the Galilean manifold to examine the weak-field 
approximation and we obtain Galilean differential equations through the 
procedures that usually describe gravitational waves in GR. This will result 
in Poisson, diffusion, and Schrodinger equations for the fluctuation field. We 
shall see in Section [3] that it is possible to accommodate MOND |10|-|13| 
within the Galilean gravity. Henceforth, we consider a 5-dimensional Rie- 
mannian spacetime; that is, a 5-dimensional differential manifold where an 
invertible symmetric metric tensor is defined. Thus we can define quantities 
usually encountered in general relativity: a symmetric connection, given by 
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the Christoffel symbols; covariant derivatives; curvature, or Riemann tensor; 
Ricci tensor, and the like. The main advantage of the Galilean covariance 
method is that these quantities have a form similar to GR, but lead to 
Galilei-invariant expressions through a natural definition of the ansatze. 

In Section [21 we consider a Lagrangian which produces equations that 
are linear in the second-order derivative of the metric. We also discuss the 
weak-field approximation with local Galilean invariance, and obtain field 
equations. We reduce the number of degrees of freedom in the perturba- 
tion field by using the trace-reverse perturbation tensor and gauge-fixing 
conditions. These results are employed in Section [3l where we examine 
three ansatze: the perturbation field does not depend on x 5 , which leads 
to a Poisson equation (Section 13. ip ; the complex ansatz, which provides a 
Schrodinger equation (Section I3.2() : and the real ansatz, which gives a dif- 
fusion equation (Section 13. 3p . The ansatz of Section [37T1 may lead to the 
AQUAdratic Lagrangian (AQUAL) field equation for MOND theory, and 
the anstze in Section 13.31 suggest a different version for the AQUAL field 
equation. We also discuss solutions of the Schrodinger equation in Section 
13.21 We make concluding remarks in Section HI 



2 Weak-field approximation in Galilean gravity 

Typically, the Lagrangian that leads to the field equations for the point- 
dependent metric is separated into two parts: Co, which depends on the 
metric, and £ m atter; which describes a matter field. The variational principle 
leads to HUES! 

$£-0 _ V~9 rpfj,u rp^v _ 2 ££matter 

5g^u 2 y/=g Sg^ 

where g = det(— g^), and T^ v is the mass- energy-momentum tensor. Hence- 
forth, greek indices, (jl, v, etc. run from 1 to 5, whereas latin indices, j, k, 
etc. denote the usual three spatial components. 

Firstly, we consider the gravitational Lagrangian to be the Einstein- 
Hilbert Lagrangian, 

A) = V~9 R, 

where R is the scalar curvature. This Lagrangian leads to Einstein's equa- 
tions [H], 

Guv = Rfiv — T.9iivR = kT^ u , (4) 
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where 

k = 8vrG, (5) 

and G is the gravitational constant. Note that we work with c = 1. 
The following relations are equivalent to Eq. 

Q -T -In T \> 

u iw flU 3 SAW ' 

where T = T^. Note that the factor 1/3, instead of the usual 1/2, follows 

from the trace of Eq. ([4]) which leads to R^ u = k (t^ v — g^uj^^T^j , with 

D = 5 instead of D = 4. 

The weak-field, or linearized- gravity, approximation is described by the 
perturbation, or fluctuation, tensor field h^, defined by 

9n» = V^u + hftv, (7) 

where | | < < 1 . Thus is the first-order correction to the Galilean met- 
ric, r}^ in Eq. ([3]). The curved spacetime described by g^ u is a perturbation 
of the flat (4 + l)-dimensional Minkowski manifold. 

Consider a finite, global (i.e. x-independent) Galilean transformation, 
in the form of Eq. ([2]), 

where kp k^n^ = r] pa , and n^ is the flat-space metric given in Eq. ([3]). 
The metric transforms as 

I Ox^ dx a Ik— WP I k —\\ a 

= A/V (r^ + V) = ^ + A/A/ V- 

We have used the property (A _1 ) p v = k u p . 

Since we have, in the new frame, g' ^ u (x') = rj^ + h' ^ (x 1 ) , we observe 
that hnv transforms like a tensor under Galilean transformations, 



U> — A PA VU 



From Eq. ([4]), we find 
1 



G 



2 



kT^, (8) 
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where the indices between parentheses are symmetrized; that is, Ar^ 
^(A^ + A VI1 ). 

We find, from Eq. ([6]), the equivalent result: 

1 



kS, 



(IV 



- \ 2 dpd^hj - d^dyh - Uh^ 



(9) 



Since R is of first order in h^ u , T^ u (or S^ u ) must also be of first order. This 
explains why we use 77^, instead of g^ u , in Eq. dH) (or Eq. ©). In this 
approximation, the interpretation of the components of the mass-energy- 
momentum tensor is manifest. As explained in Ref. [3], the components T{j 
are the density of momentum flux of the matter field; T§j is the field mo- 
mentum density; T±j is the density of energy flux; — T45 the energy density, 
and — T55 is the density of mass of the field. The component T44 does not 
have a clear physical interpretation and it seems to be simply a remnant of 
the extended manifold. 

In order to obtain a formal solution of the field equations, let us work 
with the harmonic gauge condition (see Section 7.4 of Ref. [H]): 

<rr£„ = \rr (2d v h\ - d a h) = 0. 

If we substitute this into the linearized field equation, Eq. ([9]), and then 
apply a Fourier transform, we obtain 

KSpvip) = - p p p p hp V {p). (10) 

By solving Eq. (|10p for hp U (p), and performing an inverse Fourier trans- 
form, we obtain 



h 

where 



2k J afx'G (x - x') S aT (V) , (11) 



G(x-x') = J afp e ipp{x - x ' r ^ = J afp e^-xj 



(p 2 - 2p 4 p 5 ) ' 



Note that this propagator was used in a different context in Ref. |16j . 
By integrating this expression for the Green's function, we obtain the field 
components as integrals of the energy-momentum tensor, Eq. (jlip . 

For later convenience, let us define the trace-reverse perturbation tensor: 

hnu = hp V - \rjpvh, (12) 
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so that the terms involving the trace of hn V do not appear in the field 
equation. Then we can write Eq. ([8]), i.e. the field equations without gauge 
fixing, as follows: 

KTp V = d p d {p h u) p - -rj^ d p d a h p - -Dh^. (13) 

Furthermore, with the Lorenz gauge condition, d p h u p = 0, we obtain 
the general field equations 

kT pu = ~ [v 2 - 2d±d 5 ] v. (14) 

This equation takes various forms, depending on different hypotheses con- 
cerning the functional form of the field h pu . Specific ansatze are discussed 
in Section 

(About the gauge fixing, let us comment that, on the (4 + l)-dimensional 
manifold, the rank-2 symmetric tensor h pu comprises 15 independent com- 
ponents. The Lorenz condition, d p hif = 0, introduces five more constraints, 
and the choice of coordinate system implies five additional restrictions. Thus 
the linearized tensor has only five independent components: hn, /123, /133. 
This will be transparent in Eq. (|27p of Section 13. 2\ when we discuss the so- 
lutions of the Schrodinger equation. In discussions of gravitational waves in 
five-dimensional GR, one obtains the same number of independent compo- 
nents for h pv (for instance, see Eq. (6) in Ref. |18j). Let us observe that the 
harmonic gauge condition and the Lorentz gauge condition are equivalent; 
subsequently, we choose to use one or the other at convenience.) 

3 Three ansatze for the perturbation field 

Hereafter we examine three specifc ansatze of h pu as a function of x 5 , and 
we examine Eq. (I14p . in the presence of matter-energy fields. 

3.1 Poisson equation 

If huv does not depend on x 5 , so that 

= ^(x,^ 4 ), (15) 
we see that Eq. (fT4"l) reduces to the Poisson equation: 

V 2 V(x,<) = -2kT^. (16) 
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In vacuum, T^ u — 0, it simply becomes the Laplace equation. Notice that, 
from the point of view of the Galilean approach, this ansatz means that we 
are considering hav as being a massless field. 

We already mentioned in Section [2] that — T55 is the mass density, 



T 5 5 = P- (17) 



Then, with Eq. (JTTD and 



/155 = Mn, (18) 
we find that Eq. (|16|) becomes the Newton's gravitational law: 

X7 2 cj) N = AirGp, 

where 4>n is the Newtonian potential, and k is defined in Eq. ([5]). 

Therefore, if we keep Eq. (fTT|) . but, instead of Eq. (fTHj) . we consider /155 
given by 

Vh 55 =4/i(x)V0, x = |V0|/a o (19) 



then Eq. (|16p gives the Poisson equation utilized in the AQUAL formalism 
[IHE]: 

V • [fx (x) V0] = 4^Gp. (20) 

The Milgrom's transition function [i depends only on the first derivative 
of the field (j), which is the gravitational potential that drives the motion; 
that is, a = — V</>. According to MOND, bodies subjected to gravitational 
force that move with an acceleration smaller than ao should have a dy- 
namic behaviour different from Newtonian mechanics. AQUAL reproduces 
exactly the MOND formula, /i(x)a = — V^tv, where 4>n is the Newtonian 
gravitational potential. In MOND models, [i can, for instance, consist of the 
"standard function", [i{x) = x/y/l + x 2 , with x = |a|/ao, which exhibits the 
additional fundamental constant of the theory: the modulus of the accel- 
eration, ao ~ 10~ 10 m/s 2 . For the standard /^-function given above in the 
large-^^- limit, we expand as follows: 

Va 55 ~ , 4V ^ , 
Vi+K/iv^i) 2 

= ^ f 1 - 1 (r%) 2 + 1 (r^i) 4 - ^ (r#^) 6 + • ■ 

which approaches the Newtonian limit when ao << |V0|. In this limit we 
can approximate /155 by 

^55 = 4/i (x) (j), 
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dfi 



since Vn(x) - 7E7 ^ m 
0, so that V/155 = An(x)V(j) + ^4>V n(x 



(V0 • V) V<p, and lirn^oo ^ = lim. 



3 V(i+^' 2 ) 3 

4//(x)V(/>. This identification 
cannot be done in the small-— limit, where we would obtain 

ao ' 



Vh 



55 



4V0 |V<ft|/a 

-v/l+(|V<A|/a ) : 

[Wj 
ao 



4V0 



1 /|v^l\ 3 

2 ^ a J 



_j_ 3 /IV^IV 
8 ^ a J 



_5_ 

16 



ao 



+ 



which approaches the Newtonian limit when |V0| << ao- in this case /155 
cannot be integrated and in order to obtain consistency with the AQUAL 
theory we should choose 

V/155 ~ 4// (x) Vcj) + V x f , 

where f is a function to be determined. This is in accordance to what is 
expected from the AQUAL version of MOND where the theory is determined 
up to a curl of an arbitrary vector field [13] . 

The same can be done for the "simple function", fi(x) = x/(l + x). In 
the small-^^ limit, we write V/155 as 



Vh 



55 



4V0 \V<t>\/a 
l+|V0J/a o 

4V</> "JM 



ao J ' \ ao J \ a J ~ r 



which approaches zero in the limit |V0| << ao- in this case, the /155 = n(x) 4> 
is also valid in first order of approximation. The simple ^-function, in the 
large-^^- limit, is obtained by expanding V/155 as 



Vh, 



55 



4V</> 
l+a /|V</>| 



4V0 



1 



|V0| 



+ 



ao 

|V0| 



ao 



+ 



which approaches the Newtonian limit when ao << |V^|. Here a curl of a 
vector field also should be added in order to integrate /155. 



3.2 Schrodinger equation and its wave solutions 

We now turn to the case where the gravitational perturbation is a massive 
field. Since p§ is an invariant of the Galilean algebra 116] . the natural 
ansatz is given by 

(21) 



-imx 



h/MU (x,X 4 ) . 



s 



A similar situation occurs in the Galilean analyzes of the Proca field [3]. 
In fact, if we follow our recent paper [T7], and consider both the positive 
and negative mass contributions to Eq. (|2ip . it amounts to adding the 
complex conjugate to Eq. (f2Tj) . so that h^ u is real. This is analogous to 
the treatment of plane waves discussed in Section 10.2 of [H]. The physical 
field (x, x 4 ) should be real as well. With h^ u defined as above, we see 
that Eq. (TH")) becomes 

- 2kT^ u = V% w + 2mid 4 V- ( 22 ) 

This is an inhomogeneous Schrodinger equation for each tensor component. 

Next, we find a solution for the 'free-particle' situation, or, equivalently, 
the vacuum case, T^y — 0. This leads to the Schrodinger equation, 

If we write the complex ansatz of Eq. (|2ip for h^ v in the form, 

h^ = e- imx5 T{x A )X^ v (x), 

then we attain separation of variables in the Schrodinger equation. 

If we denote the constant of separation by E, the equation for T (x 4 ) is 

d 4 T (x 4 ) +iE T (x 4 ) = 0. 

Its solution is 

T(x 4 ) =e~ iEx \ (23) 

The equation for X^ v (x) is 

V 2 X^ (x) + 2mE AV (x) = 0. (24) 
The Fourier transform of this equation is 

J d 3 pe ip -* [-p 2 + 2mE] (p) = 0, 

which, based on the independence of the plane waves for each value of p, 
gives the dispersion relation for a non-trivial solution: 



This is the non-relativistic energy of a free-particle, as expected, since we 
have chosen T^ v = 0. 

With the dispersion relation, Eq. (|25p . the solution of Eq. (|24p is given 
by an amplitude Aa V times the plane wave: 

■^■ixv (x) = e p 

The equation for X^ v (x) is solved in the particular gauge given in Section 
[2J that is, hi^ = /144 = /145 = h = 0. This leads to 

X^i = X44 = X45 = 0, 

and 

V 2 Xij (x) + 2mE Xij (x) = 0, 
V 2 X i5 (x) + 2mE X i5 (x) = 0, 
V 2 X 55 (x) + 2mE X 55 (x) = 0. 

The first equation above can be solved in the transverse-traceless ansatz, 
usual in the (3 + l)-dimensional spacetime. 

The Fourier transform of the Lorenz-like gauge condition leads to 

i Pi h i4 - LE/1 44 = 0, 
i Pi h ij - iEh 4j = 0. 

Then, we have 

?44 _ p»fe' 4 

n ~ E..> 
= 

By substituting the second equation into the first one, we obtain h 44 in 
terms of h? 3 : 

?44 _ PiPjh Jl 
~ E 2 ■ 

Consider a single plane wave; that is, consider a direction of propagation 

P 

n = - — - . 

IpI 

Furthermore take the direction of propagation along the z-axis; that is, 

Pi = P fa . 

With this restriction, some simplifications follow: 

h 44 = ^/j 33 
" E 2 
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and 



-h 3j 
E 



Together with the condition, 

h = h n + h 22 + h 33 = 0, (26) 

this gives 



TT 



Note that Eq. (|26p is also used in the standard relativistic treatment and 
it is called "traceless condition" in that relativistic context. The label TT, 
attached to huv, indicates that we have chosen the Galilean analogue of the 
transverse-traceless gauge. In terms of the amplitudes, we obtain 



/ 


hi 




hn 





fh 13 


\ 




hi2 


- (hn + fr 33 J 


h23 





1^23 






hi3 


h>23 


h33 





1^33 























V 




%h 23 


f^33 





|^33 


/ 



( 


A n 


A l2 


A 13 





1^13 


\ 




A 12 


- (A 11 + A 33 ) 


A 23 





1^23 






An 


A 23 


A 33 





^33 



























1^23 


f^33 







/ 



-i(mx 5 +Ex 4 '- 



(27) 



-pz I 



This is a wave- type solution for the Schrodinger-type equation in the transverse- 
traceless gauge. 

This suggest the following identifications: 



LO 



E, 



and k 



2tt 

T 



p, 



with phase velocity 



v 



LO _ E 

k p 

This result is also obtained in the usual quantum mechanics, described by 
the Schrbdinger equation. 

The analogue of Eq. (|27p. in the usual four-dimensional GR, exhibits 
only two degrees of freedom. We can reduce our five degrees of freedom to 
only two independent entries by setting 



A 



13 



123 



13.3 



0. 



(28) 
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in which case only An and A\2 will remain in the matrix h * . 

The meaning of the extra components in the 5-dimensional spacetime is 
still to be clarified. We remark that the additional degrees of freedom are 
associated to the massive character of the weak field solution, manifested by 
the dependence on the x 5 -coordinate. 



3.3 Diffusion equation 

Now we analyze an ansatz analogous to the previous one, but avoiding the 
complexification. The gravitational field is still supposed to be massive, 
since the dependence with the fifth-coordinate is kept. If is chosen as 
follows, 

V = e- mx V( x > a;4 )> (29) 
then Eq. (Tl4]) becomes an inhomogeneous diffusion equation, 

- 2kT^ = V 2 V + 2md 4 V, (30) 

with — ^ being the diffusion coefficient. Note that Eq. (fT9|) . substituted 
into Eq. ([30]) produces an equation similar to Eq. ([20]) but with an addi- 
tional time-derivative of /155, which is expressed in terms of (j> by integrating 
Eq. (USD. 

Notice that Eq. ([30]) is obtained when a gauge is fixed. If we suppose that 
the gauge is not fixed and that the field is static, i.e. i)^ (x, x 4 ) = rj^ (x), 
then, from Eqs. ([1~3]) and ([29]) . we find 

KT44 = — -V 2 /l44, 

kT 55 = -- (2mdihj + V 2 /i 55 ) , 
kT 45 = ^ (mdihl + V 2 /j 45 ) . 

The equations in which d{h^ and dih 5 l are non-zero constants reproduce the 
Poisson equation with extra terms analogous to the cosmological constant. 
This particular choice leads to a Newton-Hooke force; that is, a Newtonian 
potential plus a (linear) Hooke-like contribution [19] , showing that the grav- 
itational constant should be related to the mass of the gravitational field. 

Note that it is not possible to obtain such an equation in /owr-dimensional 
GR. In principle, it could be done in /iue-dimensional GR, but this would be 
far more artificial than the Galilean context utilized here. Further comments 
are in Section HI 
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4 Concluding remarks 



In this paper, we have retraced the procedures normally leading to the in- 
vestigation of gravitational waves in GR, but we have replaced the Lorentz 
metric with the Galilean metric, Eq. ([3]), so that the gravitational theory 
considered is locally Galilean invariant. The reason for doing this is that the 
Galilei group underlies the appropriate kinematics for low-energy, or non- 
relativistic, phenomena; therefore, the Galilean symmetry must be taken 
into account when non-relativistic objects are causing the gravitational per- 
turbations. 

In this paper, we examined the weak-field approximation with local 
Galilean invariance with a (4 + l)-dimensional Minkowski manifold which 
leads to Galilei-invariant equations, when the physical quantities are prop- 
erly introduced. While there exist 5-dimensional models of GR which have 
no connection with Galilean symmetry, there are a few reasons to consider 
Galilean gravity. For instance, the definition given in Eq. (119p . even though 
permitted in five-dimensional GR, is better justified in the Galilean context 
because Eq. ([T7]) associates the mass density to the component T55 of the 
energy-momentum tensor. Another feature that distinguishes 5-dimensional 
GR and the Galilean formalism employed here is that the ansatze, given in 
Eqs. (|15p . (I29p and (I2ip . appear in a natural way within the Galilean frame- 
work |16| . whereas no such motivation occurs in GR . 

The three ansatze considered in Section 3 lead to a Poisson equation, a 
diffusion equation, and a Schrodinger equation. The Poisson equation allows 
us to introduce the AQUAL field equations which describe the MOND the- 
ory. This definition and the real anzatz lead us to suggest a modified version 
of the AQUAL field equations. Whether this equation will be relevant in 
cosmology remains to be explored. 
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